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, ^ ' We present a generalized Schmidt decomposition for a pure system with any number of two-level 

\^ \ subsystems. The basis is symmetric under the permutation of the parties and is derived from the prod- 

I— I . uct state defining the injective tensor norm of the state. The largest coefficient quantifies the quantum 

pi • correlation of the state. Other coefficients have a lot of information such as the unentangled particles 

' as well as the particles whose reduced states are completely mixed. The decomposition clearly dis- 

^ ' tinguishes the states entangled in inequivalent ways and have an information on the applicability to 

O^. the teleportation and superdense coding when the given quantum state is used as a quantum channel. 

, Keywords: multipartite entanglement, separability problem, geometric measure, nonlinear eigenproblem 

> ■ 

o : 

^ . I. INTRODUCTION 

^ ■ The Schmidt decomposition for bipartite systems[l] is a very important tool in quantum information and 

Q , quantum computing theories. It shows whether two given states are related by a local unitary transformation 



or not[2], which states are applicable for perfect teleportation [3] and superdense coding[4], and whether it 
is possible to transform a given bipartite pure state to another pure state by local operations and classical 
communications [5]. Many substantial results have been obtained with the help of the Schmidt normal form 
and its generalization to the multipartite states is a task of prime importance[6, 7]. 

Carteret et al developed a method for such a generalization for pure states of an n-party system, where 
the dimensions of the individual state spaces are finite but otherwise arbitrary [8]. The idea of this method 
is the following. First, one finds the product vector that gives maximal overlap with a given state func- 
tion. Next, one considers product vectors whose constituents are orthogonal to the constituents of the first 
product vector and finds among them the product vector that gives maximal overlap with the state function. 
Continuing in this way, one finds a set of orthonormal product states and presents the state function as a lin- 
ear combination of these product vectors. Since the first product vector is a stationarity point, the resulting 
canonical form contains a minimal set of state parameters. 
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Surprisingly, an arbitrary multipartite quantum state has many canonical forms. The reason is that 
the stationarity equations(SEQ) defining stationarity points are nonlinear equations and in general have 
many solutions of different types [9]. For example, in the case of generic three qubit states there are six 
product vectors that are solutions of SEQ[10]. Moreover, some highly entangled states are surrounded by 
the one-parametric set of equally distant product states that have the same maximal overlap with the state 
function[ll]. All of these product vectors are equally good in a sense that all resulting canonical forms use 
minimal number of product vectors from a factorizable orthonormal basis. Therefore, it is necessary to make 
a physically motivated choice of the product vector that yields the best factorizable basis. More precisely, 
the problem is the following. Each of bases suggests its own set of invariants which are coefficients of 
the expansion of the state function. And the task is to specify the set of invariants that gives an effective 
description of the quantum state. 

Our approach is based on the injective tensor norm of quantum states[12]. The overlap of a given state 
with any complete product state attains maximums at several points. These maximums are eigenvalues and 
corresponding product states are eigenvectors of SEQ. The injective tensor norm takes the value of one of 
these maximums and thus makes a unique choice of the eigenvector. This choice classifies pure quantum 
states by types as follows. Consider three qubit pure states. There are six eigenvectors and eigenvalues of 
SEQ and therefore the injective tensor norm can have six different expressions. Since it is a single valued 
function, each of expressions should has its own range of definition in which they are deemed applicable. 
Consequently, each of solutions of SEQ has its own applicable domain depending on state parameters in 
which they give the relevant eigenvalue. As a result, the space of the state parameters is split into six regions 
and each region has its own specific eigenvector defining the injective tensor norm. Each of three qubit pure 
states should belong to the one of six regions of the parameter space except states that lie on joint surfaces 
separating different regions. This classification provides a geometric picture of entangled regions. 

The main objective of this work is to construct the generaUzed Schmidt decomposition(GSD) for generic 
three qubit states. We analyze the complete set of SEQ solutions, specify the relevant solution for a given 
state and calculate GSD coefficients expUcitly in terms of the state parameters. In this way we obtain a 
new set of invariants of three qubit states. Since there are different sets of local invariants[6, 13], one may 
wonder what is the advantage of the obtained coefficients. For this purpose we study the properties of the 
constructed GSD. 

The expansion coefficients exhibit the physically significant properties of pure states. The largest coef- 
ficient g is the injective tensor norm of the state. It is a very useful quantity and defines some entanglement 
measures [9, 14, 15, 16, 17]. Apart from that, we will present a rigorous proof that if a three qubit pure 
state has a completely mixed one qubit reduced state, then = 1/2. This sheds some new light on the 
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applicability of the quantum state to the teleportation and superdense coding. There is a conjecture that 
the criterion for perfect quantum teleportation is = 1/2[18]. Now this conjecture can be reformulated 
as follows: a pure state can be used for perfect teleportation and superdense coding if it has a completely 
mixed state as a reduced state. The other expansion coefficient, say h, has an information on the presence 
or absence of an unentangled particle in a given quantum state. We will show in the following that /i = is 
a separabihty criterion for pure states of a general multi-qubit system [19]. 

This paper is organized as follows. In Sec. 11 we specify GSD for multiqubit systems. In Sec. Ill we 
analyze general properties of GSD in the case of three qubit systems. In Sec. IV we construct GSD for 
three qubit W-type states. In Sec. V we construct GSD for a class of multiqubit W-type states. In Sec. VI 
we construct GSD for Greenberger-Home-Zeilinger(GHZ)-type states. In Sec. Vn we make concluding 
remarks. 

II. GENERALIZED SCHMIDT DECOMPOSITION 

Consider n-partite pure systems with the Hilbert space H = Hi (8) H2 (8> • • • (8) Tin- The injective tensor 
norm g{ip) of a given n-partite pure state is defined as 

=SUp|(xiX2---Xn|V')|, (1) 

where the supremum is over all tuples of vectors \xk) € Hk with \\xk\\ = 1 [12]. The nearest product state 
\q) = \qi Q2 ■ ■ ■ Qn) must satisfy stationarity equations [9, 20] 

{qiq2---qk---qn\'<P) = gWk), k = l,2,---n (2) 

where the caret means exclusion. These equations can be regarded as nonlinear eigenvalue problem: all 
eigenvalues g and eigenvectors \q) satisfy it for a given state function In general, Eq.(2) has several 
different solutions. Hereafter we consider only the solution corresponding to the eigenvalue of largest value. 
This solution, known as dominant eigenvector, defines the injective tensor norm of a given quantum state. 

Consider now n-qubit system. For each single-qubit state \qk) there is, up to arbitrary phase, a unique 
single-qubit state \pk) orthogonal to it. From these single-qubit states \qk) and \pk) one can form a set of 2^* 
n-qubit product states which form a basis in the full Hilbert space H. Any vector £ H can be written 
as a linear combination of vectors in the set. Then from SEQ (2) it follows that all the coefficients of the 
product states \qi - ■ ■ qk-iPkQk+i • " ■ Qn){k = 1, 2, • • • n) are zero. Thus any pure state can be written in 
terms of 2" — n product states. Furthermore, the phases of vectors \pk) are free and we can choose them so 
that all the coefficients tk of vectors \pi ■ ■ -Pk-iQkPk+i ■ ■ ■Pn){k = 1, 2, • • • n) be positive. Still we have 
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a freedom to make a phase shift \pk) —>■ e^'^^/^^^^^lpfe) which remains unchanged and g. We use this 
freedom to vary the phase ip of the component e^'^h\p\p2 ■ ■ ■ Pn) > is understood) within the interval 
—it/ (n — 1) < < 7r/(n — 1). 

Thus the decomposition has n + 1 real and 2" — 2n — 1 complex parameters. After taking into account the 
normalization condition, one can show that 2"+^ — 3n — 2 real numbers parameterize the sets of inequivalent 
pure states [21]. 

Theorem 1. The A;th qubit is completely unentangled if and only if h{^) = and ti(V') = for i 7^ A;. 
Proof. Suppose first qubit is completely unentangled and its state vector is \q\). We have |^) = <8) 
l'^'). Let the product state \q2qz ■ ■ ■ Qn) be the nearest state of Then GSD of l'^') takes the form 

n 

IV'') = 9'\q2q3 • • • 9n) + ^ |P2 • • • Pi-mPi+l ■■■Pn) -\ h e'"^' h'\p2P3 ■■■Pn)- (3) 

i=2 

Since the nearest state of the state is, up to a phase, the product state \qiq2 - ■ ■ Qn)^ then g{ijj) = 
g', h{ijj) = 0, ti(V') = h' and ti = 0,i = 2, 3...n. The inverse is also true. From h{ijj) = and ti{ip) = 
for i 7^ 1 it follows that all the terms in GSD which do not contain \qi) vanish and 1-0) = l^i) <8) 
Similarly, theorem is true if any other qubit is unentangled. 

Consider now n = 2 and n = 3 cases. For simplicity we will use notations |0j) and for vectors 
\qi) and \pi) respectively. Also we will omit sub-indices i whenever it does not create misunderstanding. In 
the case of two qubit states the expansion reduces to the Schmidt decomposition lip) = g\00) + h\ll) with 
g > h>0. Consider three-qubit case. 

m. THREE QUBIT STATES. 

Decomposition takes the form 

IV-) = g\000) + ii|011) + i2|101) + talllO) + e^^/i|lll). (4) 
The coefficients should satisfy conditions 

g>ms^{ti,t2,t3,h), ti>0, t2 > 0, ts > 0, h > 0, ~ < ip < (5) 

These conditions do not specify GSD uniquely. Eq.(4) is the GSD normal form of the state |^) if and only 
if g is the injective norm of the state 1-0). There are highly entangled states which can be written in a form 
of Eq.(4) in six different bases. One basis, where the largest coefficient is injective tensor norm of the state 
and corresponding product state is the dominant eigenvector of SEQ, gives true GSD. The others, where 
the largest coefficient is not injective tensor norm and corresponding product states are eigenvectors but are 
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not dominant eigenvectors, do not. The example with W-type states, which is given in the next section, 
illustrates this more clearly. 

Now we formulate a theorem which shows whether a given three qubit pure state has a completely mixed 
one qubit reduced state or not. The proof of Theorem 2 uses a particular lower bound on g given in the next 
section. In this reason we present the proof in Appendix. 

Theorem 2. One-particle reduced state of kHa qubit is completely mixed if and only if 

tfc = and = 1/2. 
Consider now several interesting examples. 

IV. W-TYPE STATES 

Our first example that we shall discuss in detail is a family of four-parametric W-type states [22] 

IV-) = a|100) + 6|010) + c|001) + d\l\l). (6) 
Without loss of generahty we consider only the case of positive parameters a, 6, c, d. 



A. Solutions of stationarity equations 

Stationarity equations (2) have six different solutions for each of these states [10]. Four solutions are 
simple and represent four terms on the right-hand side of Eq.(6) 

|gig2g3) = |100), g = a; (7a) 



1919253) = |010), g = h- (7b) 
\q1q2qz) = |001), g = c\ (7c) 



|gi9293) = g = d. (7d) 

These four solutions owe their existence to the fact that the states Eq.(6) have a hidden symmetry. Namely, 
there exist local unitary transformations that interchange the positions of any of two coefficients in Eq.(6). 
Due to this symmetry any entanglement measure is invariant under the permutations of all of the state 
parameters a, 6, c, d. 
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The fifth nontrivial solution is 
where 

r„ = a(6^ + c'^ + (f- a^) + 26cd, = b{a^ + c'^ + (f- h^) + 2acd, rc = c{a^ + h'^ + d^- c^) + 2ahd, 



Td = d{a^ + 6^ + - d^) + 2a6c, L = y/{ab + cd) (ac + 60?) (ad + be) 

and S is the area of the cychc quadrangle with sides (a, b, c, d). Note that the quadrangle is not unique as 
the sides can be arranged in different orders. But all these quadrangles have the same area given by Heron's 



formula S = ^/{s — a){s — b){s — c)(s — d), where s = {a + b + c + d)/2 is semiperimeter. 
The sixth solution is obtained from the solution (8) by the interchange 

d -d, rd -rd- (9) 

Due to symmetry one can make one of the following interchanges instead 

a -a, ra -ra, 
b -b, rb -rb] 
c -c, Tc -rc 

and obtain the same solution. However, it does not mean that any of parameters becomes negative. All the 
parameters are positive and the interchange (9) shows the existence of a pair of dual solutions. 



B. Dominant eigenvector 

We have six different product states \q1q2Q3) and therefore we can form six sets of factorizable bases 
from vectors \qi) , \ q2) , ll^) and \pi), \p2), Ips)- Since all these product states satisfy SEQ, all the expansions 
of the state function in these bases give the canonical form Eq.(4). Therefore there are six factorizable 
bases where the state function Eq.(6) has a canonical form Eq.(4). The question at issue is whether there 
are effective expansions or not. If it is indeed the case, then what are the effective expansions which have 
physically meaningful coefficients? 

The injective tensor norm g has six different expressions. Since it is a single valued function, each 
expression has its own applicable domain depending on state parameters and these appUcable domains are 
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split up by separating surfaces. Consequentiy each of solutions (7) and (8) has its own appUcable domain 
depending on state parameters in which they are dominant eigenvectors. Thus four-dimensional sphere 
given by normalization condition + P + + cf = 1 is split into six parts and each part has its own 
dominant eigenvector 

The above classification suggests a natural choice of a factorizable basis. We choose the product state 
\Q1Q2Q3) that is dominant eigenvector in a given part of sphere. Next we form product states from con- 
stituents \qi) and orthogonal to them vectors \pi). Thus it is necessary to separate the validity domains and 
to make clear which of the solutions should be applied for a given state. It is a nontrivial task and for a 
thorough analysis we refer to Ref.[10]. Here we present some necessary results. 

In highly entangled region parameters (a, b, c, d) form a cychc quadrilateral, injective tensor norm is 
expressed in terms of the the circumradius of the quadrangle and corresponding solution is Eq.(8). In 
slightly entangled region injective tensor norm is the largest coefficient and corresponding solutions are 
Eq.(7). Also there are states in between for which both formulae are valid. These states, called second 
type shared quantum states, separate slightly and highly entangled states and can be ascribed to both types. 
Another specific states, called first type shared quantum states, are those for which injective tensor norm is 
a constant and is defined hy = 1/2. These states allow perfect quantum teleportation and superdense 
coding scenario. 

C. Highly entangled W-type states 

Highly entangled region is defined by inequahties 

ra>0, n>0, rc>0, rd>0. (10) 

The constituents of the dominant eigenvector are given by Eq.(8) and the task is to present the state function 
as a linear combination of product vectors 

IV') = 9\QiQ2q3) + ti\qiP2P3) + ^21^152^3) + ^31^1^293) + e''^h\pip2P3), (H) 

where {pi\qi) = 0, i = 1,2,3. 

The calculation of the coefficients requires advanced mathematical technique. One has to factorize 
polynomials of degree ten. We would like to suggest a simple way. First one convinces oneself that each 
factor below is a root for the corresponding polynomial and next finds the proportionality coefficient in 
some particular case. The derivation of h is the most complicated out of all coefficients and one can use the 
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hint: ifa = b + c + d, then Vb = rc = = —ra- The resulting answer is 

- Jl f - -^'-i , _ Lr2 _ Lrs _ ^ i, _ VfanrcTd 

^~ 2S' ^~ 4S{ad + bc)' ^~ AS{bd + ac)' ^~ 4S{cd + ab)' 2' ALS ' ^ ^ 

where ri, r2, are the lengths of the Bloch vectors of first, second and third qubits 

= + - - d^l, r2 = la^ + - 6^ _ rfS,^ ^ |^2 + ^,2 _ ^2 _ ^2| (^3) 

In fact, this set gives a fruitful description of the state. The invariant g is expressed in terms of the circum- 
radius of the cychc quadrangle a, b, c, d and gives geometric and Groverian entanglement measures of the 
state. First type shared states are defined by rir2r3 = which is to say that the reduced state of one of 
particles is completely mixed. It is quite obvious that one of coefficients U must vanish for these states. 
On the other hand if = 0, then = 1/2 and the corresponding state allows teleportation(and dense 
coding) scenario. For perfect teleportation the receiver should choose A;th particle at initial stage in order to 
perform the task. Thus the coefficients ti contain an information on the applicability to the teleportation and 
precisely indicate which particle the receiver should choose. Second type shared states lie on the separating 
surface and this surface is defined by raVhrcrd = 0, i.e /i = 0. We conclude that /i > for highly entangled 
states and /i = for second type shared states. 

To complete the analysis let us consider the remaining slightly entangled cases. 



D. Slightly entangled W-type states 

At least one of quantities ra, r^, Vc and should be negative in this case. Note that two of them can not 
be negative simultaneously. Indeed, 

bra + arb = 2{ac + bd){bc + ad)>Q (14) 

and therefore Va and ri, can not be negative together. Similarly, any pair of quantities ra,rh,rc, can not 
be negative simultaneously. Thus one of them should be negative and others should be positive and there 
are four possibilities. Consider these cases in turn. 
a. Solution (7a). The first case is 

Ta < 0. (15) 

The dominant eigenvector is given by (7a). In order to obtain GSD one has to simply relabel the basis of 
the first qubit |0i) ^ Then the final GSD coefficients are 

g = a, h = 0, ti = d, t2 = c, ts = b. (16) 



b. Solution (7b). The second case is 

n < 0. (17) 

The dominant eigenvector is given by (7b). Now one has to relabel the basis of the second qubit IO2) 1 12)- 
GSD coefficients are 

g = b, h = 0, ti = c, t2 = d, ts = a. (18) 

c. Solution (7c). The third case is 

To < 0. (19) 

The dominant eigenvector is given by (7c) and one has to relabel the basis of the third qubit IO3) II3). 
GSD coefficients are 

g = c, h = 0, ti = b, t-2 = a, = d. (20) 

d. Solution (7d). The fourth case is 

Td < 0. (21) 

The dominant eigenvector is given by (7d) and one has to relabel the bases of all qubits |0i) <-> | li) , IO2) 
II2), IO3) ^ II3). Resulting GSD coefficients are 

g = d, h = 0, ti = a, t2 = b, ts = c. (22) 

All of four cases can be summarized as follows. If the state Eq.(6) is sUghtly entangled, then 

i) its last coefficient ^ = vanishes and g takes the value of the largest coefficient 

h = 0, g = max(a, b, c, d) (23) 

ii) GSD takes the form 

\^) = 5|000) + tilOll) + 12|101) + 13|110), (24) 

iii) All the inequaUties (15), (17), (19) and (21) become the same lower bound on g 

g'>tl + tl + tl + 2^-^. (25) 

The obvious conclusion is that h ^ only for the highly entangled states and identically vanishes for 
the sUghtly entangled states. 
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Note that four slightly entangled cases (15), (17), (19) and (21) together with the highly entangled case 
(10) cover the fourdimensional sphere given by the normalization condition. Thus the sixth solution does 
not create its own factorizable basis. However, it does not mean at all that the solution is unphysical. The 
state Eq.(6) is a four parametric state whiles generic three qubit states should have six parameters. Thus the 
state Eq.(6) is not generic and the applicable domain of the sixth's solution has been shrunk to the separating 
surface rir2r^ = in this case. It is easy to check that the sixth solution results the same GSD coefficients 
Eq.(12) for first type shared states. But it will create its own GSD for generic states. 

Let's now examine multiqubit W-states. 



V. MULTIQUBIT W-TYPE STATES 



In this section we consider one-parametric n-qubit W-states 



a (|100 • • • 0) + |0100 • • • 0) • • • |00 • • • 010)) + 6|00 • • • 01). 



(26) 



Slightly entangled region is given by rn = {n — l)a? — < [23]. In this region the last product state 
|0 • • • 01) is the nearest separable state and 



g = b, h = 0. 



(27) 



In highly entangled region r„ > and, consequently, Sn = {n — l)^a^ — 6^ > 0. The constituent states for 
the closest separable states are respectively 

aV(n-l)(n-2)|0) + v^|l) _ ^(n - l)r„|0) + 6V^r^|l) 



ki) = --- = 



In-l 



) = 



Straightforward calculation gives 



g = il- 62)7+1/2 
where 



n-2 



, tn = A/(n - 2)r 



h = bV n — 1 



n-V 



(28) 



(29) 



n-2 



Expressions (29) have the same meanings as in the three-qubit case. First, r„ = forces (/^ = 1 /2. 
Second, g'^ > 1/2 and h = means the state is slightly entangled. Third, g"^ < 1/2 and h = means 6 = 
and, therefore, the last qubit is unentangled. Fourth, we conjecture that: all the states with r„ = allow the 
teleportation scenario and the receiver should choose nth qubit. In summary, suggested GSD indicates the 
applicability to the teleportation and distinguishes the unentangled particles as well as completely entangled 
particles. 
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VI. GHZ-TYPE STATES. 



Consider now the extended GHZ state [24] 

\'ip) = a|000) + 6|001) + c|110) + (30) 
Stationarity equations have two solutions[ll] 

ki) = |0), k2) = |0), |g3) = "'?^^\ 5= Va2 + 62; (31a) 

M = |l), k2) = |l), |g3) = ''°^^\ 5=v/^^- (31b) 

V C"' + M 

The state function takes the canonical form Eq.(4) in two bases. In both cases nonzero coefficients of the 
decomposition can be presented as follows 



g = u.^(V^fl^\V^^), t, = ^^, h=^-^^. (32) 



This set of GSD coefficients describes the extended GHZ-type states almost in the same way as bipartite 
systems. Since g'^ > 1/2, there is no highly entangled region for GHZ-type states. In this sense W-state 
is more entangled than GHZ-state. When the extended GHZ-state is most entangled, i.e. g"^ = 1/2, it is 
applicable for both teleportation and dense coding and the situation is same in the case of bipartite systems. 
In contrast to W-type case, there is no region where h is identically zero. Only on condition ad = be the 
canonical coordinate h vanishes. Thus if h vanishes, then the state is biseparable and again the same is 
true for two-qubit systems. The only difference from two-qubit case is that there is an extra term with the 
coefficient ts. It shows that the third particle is unentangled when h = 0. 



Vn. SUMMARY 



We have used a general way of constructing generalized Schmidt decomposition proposed in Ref.[8] and 
specified a unique decomposition for arbitrary composite systems consisting of two-level subsystems. This 
decomposition suggests a new set of local invariants. Three qubit pure states have following six invariants 

9, h, t2, ts, h, ip. 

We have considered general three-qubit and W-type n-qubit systems whose injective tensor norms were 
already derived analytically. All the invariants have been calculated explicitly and expressed in terms of state 
parameters a,b,c, d. It is shown that they provide a profound information on the quantum states. The largest 
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coefficient g gives two entanglement measures and together with the last coefficient h clearly distinguishes 
the states entangled in inequivalent ways. Namely, for W-type states there is a region where the function 
/i(a, 6, c, d) is identically zero and there is no such region for GHZ-type states. Furthermore, isolated zeros 
of the function h{a, b, c, d) indicate the appearance of the unentangled particles. The coefficients U are 
related to Bloch vectors and reveal the existence of completely mixed reduced states. Thereby they show 
whether or not a given state is appUcable for perfect teleportation(and dense coding) and precisely indicate 
which particle the receiver should choose at initial stage in order to perform the task. 

We have derived a lower bound on injective tensor norm Eq.(25) for slightly entangled states(h=0). A 
natural question now arises: is it possible to derive such a lower bound for arbitrary three qubit states? This 
question is closely related to the following problem. For every measure, including geometric, there must be 
a set of states which are maximally entangled, or at least approaching some lower bound. And the problem 
is to find these states explicitly. If one derives a strong lower bound on g, then one can see when the lower 
bound is saturated and obtain maximally entangled states. 

A lower bound on g for generic three qubit states can be derived by making use of GSD form (4). Indeed, 
consider a state \iprest) = ti|011) +t2|101) + i3|110) +e*^/i|lll). Suppose \pr) is its nearest product state 
and fi{ti,t2, t3,h) is its injective tensor norm. Obviously 

g > \{pr\tp)\ = \g{pr\000) + 

This is a generic lower bound and its derivation reduces to the derivation of \pr) . Unfortunately, known 
methods do not allow us to solve SEQ for IV'rest) and find \pr). Since the phases can be absorbed in 
the definitions of basis vectors, IV'rest) is a four parametric state. It is interesting to note that the states 
of type IV'rest) are the only remaining four parametric states whose injective norm has not been obtained 
analytically so far. There is a good reason to calculate this injective tensor norm and obtain a lower bound 
on it for three qubit states. 
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APPENDIX A: PROOF OF THEOREM 2 

Proof. Suppose = and = 1/2. Equation (4) allows to express the Bloch vectors of qubits in 
terms of GSD coefficients. Denote by ri the norm of the Bloch vector of the first particle. Straightforward 
calculation yields 

rl = AhHl + {g'' + tl-tl-tl-h^f . (Al) 

Then ti = and = 1/2 together with the normaUzation condition g'^ + tf + + + h"^ = 1 force 
ri = 0. Therefore the reduced state of the first particle is completely mixed. 

Suppose now that the reduced state of the first particle is a completely mixed state. Then ri should vanish 
which is to say that 

hti = 0, g"^ + tl-tl-tl- h'^ = 0. (A2) 

Either ti = or h = 0. If h = 0, then the state Eq.(4) can be transformed to the state Eq.(6) by a local 
unitary transformation |0i) <-> In this reason both states have the same GSD and Eq.(4) coincides with 
Eq.(24). Thus the inequality Eq.(25) holds for both states. On the other hand h = forces c?^ = + *3 ~ *i 
and this is contradicted by inequality Eq.(25). Consequently ii = and the the normaUzation condition 
forces g^ = 1/2. This ends the proof of the theorem. 



